To appear in The Astrophysical Journal 



Non-Equilibrium Photodissociation Regions: 
lonization-Dissociation Fronts 

Frank Bertoldi 

Max-Planck-Institut fiir Extraterrestrische Physik, D-85740 Garching, Germany 

fkb@mpe-garching.mpg.de 

and 

B.T.Draine 

Princeton University Observatory, Peyton Hall, Princeton, NJ 08544, USA 

draine@astro.princeton.edu 

ABSTRACT 

We discuss the theory of coupled ionization-dissociation fronts produced 
when molecular clouds are exposed to A < lllOA radiation from hot stars. A 
steady, composite structure is developed, which generally includes an ionized 
outflow away from the cloud, an ionization front, a layer of photodissociated 
gas, a photodissociation front, and a shock wave preceding the photodissociation 
front. We show that the properties of the structure are determined by two 
dimensionless parameters, -0 and S, and by the Alfven speed in the preshock 
gas. For a broad range of parameters of interest, the ionization front and the 
hydrogen photodissociation front do not separate, the H2 photodissociation and 
photoionization take place together, and a classical hydrogen "photodissociation 
region" (PDR) does not exist. We also show that even when a distinct 
photodissociation region exists, in many cases the dissociation front propagates 
too rapidly for the usual stationary models of PDRs to be applicable. We 
discuss several famous PDRs, e.g., in M17 and Orion and conclude that they 
cannot be described by equilibrium PDR models. 

Subject headings: ISM: molecules - molecular processes - galaxies: ISM - 
nebulae: H II regions, Orion Nebula - shock waves 
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1. Introduction 

A photodissociation region (PDR) is tlie transition layer at the surface of a molecular 
cloud between the dense, cold molecular gas and the tenuous, warm, ionized gas; it is 
the region where the far-ultraviolet (FUV) photons in the range 912 - 1110 A are either 
absorbed by H2 (which in ~ 15% of the cases leads to dissociation) or by dust. Over the last 
decade observations have shown that a significant fraction of the energy output from massive 
stars is absorbed and reprocessed in PDRs, and they have received increased theoretical 
attention [see the review by Hollenbach & Tielens (1994), and references therein]. The 
quantitative comparison of PDR models with observations can yield important information 
on the chemical and topological structure of molecular clouds and the effects of massive 
star formation on them. 

The sharp transition where molecular hydrogen turns mostly atomic we call the 
photodissociation front (DF); it will propagate into the molecular cloud when the incident 
FUV flux increases or the self-shielding neutral gas and dust layer (i.e., the PDR) 
evaporates off the surface of the cloud. This occurs when the atomic gas layer is eroded 
by a propagating photoionization front (IF) at which the neutral gas is ionized by the 
Lyman continuum (Lyc: A < 912A) radiation and, because of its sudden heating to 10^ K, 
accelerates away from the cloud surface. Note that the most prominent PDRs, i.e. those 
receiving a strong UV flux, are usually also subject to erosion by Lyc radiation, and that 
therefore the structure and dynamics of such PDRs cannot be considered without including 
the effects of the ionizing radiation. 

Semantic differences exist in the definition of what constitutes a PDR. The acronym 
is also used for "photon-dominated region," which typically denotes a surface layer that 
reaches deeper into a molecular cloud, and in which, despite its exponential attenuation, the 
FUV radiation dominates the chemical and thermal structure of the gas. In fact, most of 
the cold ISM is dominated by FUV ionization (e.g. McKee 1989), and only a small fraction 
of molecular gas is sufficiently shielded that cosmic rays are the main source of ionization. 

To distinguish the cold, low-ionization gas from the infrared-bright surface layers 
where most of the UV radiation is absorbed, for us "PDR" shall denote only the neutral 
H layer surrounding a cold molecular cloud where H2 has recently been dissociated by 
FUV photons. The PDR is bounded by a H ionization and a H2 dissociation front. The 
photodissociation of molecules such as CO can occur upstream of the H2 dissociation front, 
and some of the characteristic emission associated with PDRs (from e.g. O, C, or C"*") may 
arise there. 



Most theoretical PDR models so far assumed the gas to be at rest with respect to 
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the DF, and most (e.g. Tielens & Hollenbach 1985; Burton et al. 1990; etc.) assumed a 
steady-state balance between H2 destruction (mainly by UV pumping) and H2 formation 
(mainly by grain catalysis). 

Time-dependent PDRs have previously been discussed by Hill & Hollenbach (1978) and 
Roger & Dewdney (1992), who studied the interaction of a spherical ionization-shock front 
with the DF in an expanding H II region. London (1978) and Goldshmidt & Sternberg 
(1995) follow the propagation of a DF into a planar cloud when the radiation source first 
turns on. They find that advancing DFs produce larger column densities of vibrationally 
excited H2 than stationary ones because grains initially absorb a smaller fraction of the 
UV photons. Wagenblast & Hartquist (1988) also studied the changing H2 level population 
when the DF first penetrates a translucent cloud. 

Giant molecular clouds are known to be clumpy, with most of the molecular mass 
confined to clouds (or clumps) that fill only a small fraction of the volume of a CMC (e.g. 
Bertoldi & McKee 1992; Wilhams & Blitz 1995). When an O star creates an H II region 
within a clumpy GMC, most of the ionization and dissociation therefore occurs on the 
surfaces of the small, dense molecular clouds. None of the previous works considered the 
interaction of an ionization-shock front with a dissociation front while both fronts penetrate 
such a spheroidal molecular cloud. In this paper we show that the flow of the gas through 
the DF, PDR, IF, and off the cloud surface can have serious effects on the PDR structure 
not only when a cloud is first exposed to the radiation, but throughout the lifetime of the 
massive star or cloud. The effects of advection on the structure and emission characteristics 
of PDRs must therefore receive special attention in future PDR models. 

In section 2 we broadly specify for which astrophysical objects we expect the IF to have 
little effect on the structure of the PDR. In section 3 we introduce the two main parameters 
that define the impact of the ionizing radiation on a cloud; the photoevaporation parameter, 
-0 is a measure for the Lyc opacity of a steady photoevaporation flow off the surface of 
the cloud, and the Str0mgren parameter, S, is a measure for the dynamical impact of the 
photoevaporation on the cloud. In section 4 we compute the propagation speeds of the 
ionization and H2 dissociation fronts and determine in which cases both fronts form a 
merged structure, so that an extended neutral HI layer can not develop. We flnd that this 
condition mainly depends on the photoevaporation parameter of the cloud and the ratio 
of the Lyc and FUV emission of the star. Section 5 compares the flow time through a 
propagating DF with the timescale within which equilibrium H2 abundances are established. 
We can thereby show in which cases transient photochemistry has serious effects on the 
structure of the PDR. In section 6 we discuss the implications of our flndings for several 
well-known objects, including the PDR behind the Orion Trapezium stars, and section 7 
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summarizes our findings and discusses their implications. 

2. Cases Unaffected by Ionizing Radiation 

The structure of PDRs is most affected by ionizing radiation where the latter is strong 
compared to the FUV. We focus our discussion therefore on molecular clouds that are 
irradiated by one or more nearby O stars. We thereby exclude objects such as reflection 
nebulae, H I regions dissociated by the diffuse FUV background radiation, and the edges of 
certain H II regions. 

2.1. Reflection Nebulae 

B stars have very weak Lyc emission but can form prominent shells of warm, dissociated 
H I around their small H II regions; NGC 2023 is a famous example of such a "reflection 
nebula" (e.g. Gatley et al. 1986; Draine & Bertoldi 1996). The dynamics of such PDRs are 
not strongly affected by ionizing radiation. 

2.2. Clumps in Molecular Clouds 

Far infrared observations of e.g. C"*" (Stutzki et al. 1988) show that the diffuse 
FUV background radiation can apparently propagate deep into clumpy molecular clouds, 
dissociate the tenuous interclump gas and irradiate the surfaces of dense molecular clumps. 
The PDRs around such clumps are likely in pressure equilibrium with the interclump gas 
and are not affected by ionizing radiation. However, if the FUV in such a clumpy cloud is 
enhanced by nearby massive stars (as in the Rosette Molecular Cloud: Cox et al. 1990) and 
is increasing due to the approach of an ionization front, time-dependent effects may become 
important; the PDRs may begin to propagate into the clumps in response to the increasing 
FUV flux. 



2.3. Edges of Spherical H II Regions 



At the edge of an H II region that expands into a homogeneous molecular cloud (i.e., an 
ideal Str0mgren sphere), the ionizing radiation is weak and the IF propagates only slowly 
into the PDR that forms the inner boundary to the dense neutral shell piled up during 
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the expansion. If photoevaporating clumps inside the H II region contribute significantly 
to its mass input (McKee et al. 1984), the edge of the H II region can even become a 
recombination front (Bertoldi 1989a). In either case the ionizing radiation does not have a 
strong effect on the structure of the PDR, and quasi-equilibrium models (Hill & HoUenbach 
1978; Roger & Dewdney 1992) are appropriate to describe the evolution of such PDRs, 
except during the initial phases of their expansion when the IF is R-type. 



3. Photoevaporation Flow Parameters: tp and S 

We now consider PDRs on the surfaces of molecular clouds which are illuminated by 
one or more O stars. We assume that the surface is convex, i.e., curved away from the star, 
so that the newly ionized gas can freely expand into the H II region (Fig.l). Let rg be the 
radius of curvature of the surface, and R be the distance from the surface to the star. The 
structure of such a photoevaporation flow, the ionization front, and the compressive effect 
on the cloud were discussed in detail by Bertoldi (1989b) and Bertoldi &: McKee (1990) for 
the case re- They showed that the photoevaporation of a spheroidal cloud is essentially 
characterized by two parameters. The dimensionless "photoevaporation parameter" -0 is a 
measure for the opacity of the evaporation flow to Lyc photons: 

^ = A m2 = 5-15 X 10^ — ^ , 1 

inR^c^ R^^ 

where ct^ ~ 2.59 x 10^^'^cm'^s^^ is the case B hydrogen recombination coefficient at 10'' K, 
C2 — 11.4 km s~^ is the typical isothermal sound speed in the ionized gas, = Tpcpc, 
R = RpcPC, and Sly = 10''^ ^'49 s^^ is the emission rate of Lyc photons by the star. Typical 
photoevaporation flows of interest have 0^1. 

Here we consider the more general case, including the limit of a planar cloud, R <^ rg. 
Let r be the distance from the ionizing source. Along the symmetry axis, the ionized gas 
streaming off the cloud surface has a density profile approximated by 

n{x) ^ n2{x/r,)-^-^ , (2) 

where x = R — r + rc is the distance from the center of curvature, and n2 = Fiy/v2, where 
Fly is the ionizing flux arriving at the cloud surface, and V2 ~ I.6C2 is the characteristic 
speed of the evaporation flow. The exponent 2.5 and coefficient 1.6 were chosen to fit the 
high-'i/) density profiles computed by Bertoldi (1989b). 

In ionization equilibrium, the radiation transfer equation along a radial line from the 
star is id 

dr l^'^^^'^^] ^ -aBn{r)'^ - n{r)aiy F{r), (3) 
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where F{r) is the Lyc flux at some distance r from the star, and aiy is the dust absorption 
cross section per H nucleus for the ionizing photons. We adopt the "on-the-spot" 
approximation, assuming that recombination to tlie ground state and the resulting Lyman 
continuum photons can together be omitted. Similarly, scattering of Lyman continuum 
photons by dust is neglected, as these scattered photons remain available for ionizing 
hydrogen. 



Interstellar dust properties vary from region to region [see the review by Mathis 
(1990)]. The extinction longward of / = 9000A appears to be "universal" (Cardelli, Clayton 
& Mathis 1989), with Ai/Nu ^ 2.60 x lO^^^cm^ (Draine 1989). For A < 9000A, the 
extinction law varies, depending on the value of Ry = Ay/ E{B — V). Diffuse cloud dust 
is characterized by Ry ~ 3.1, and the extinction cross section per H nucleus is estimated 
to be - 2.6 X 10-21 cm^ at lOOOA (Cardelli, Clayton & Mathis 1989). For dense clouds, 
however, the dust is characterized by larger values of Ry ~ 4 — 6 (Mathis 1990), and 
reduced FUV extinction. For the numerical examples in most of this paper we will consider 
dust characterized by an average dense cloud Ry 5, close to the high value of 5.5 that is 
found in the Orion (M42) nebula (§6.3). 

Grain models predict that the absorption will rise to a maximum near ~ 750A, beyond 
which it will decline (Laor & Draine 1993). Since we are mostly concerned with hot stars 
in dense H II regions, the ionizing photons near the ionization front should have average 
photon energies well above 13.6eV. Thus the effective extinction cross section for these 
photons is smaller than the maximum value at 800A. In lieu of reliable extinction models, 
we will assume that the average extinction cross section is close to that extrapolated to the 
Lyman edge, ~ 1.22 x 10^^^ cm^/H; with an estimated albedo of 0.5, the Lyc absorption 
cross section is then aiy ~ 0.61 x lO^^^cm^. 



3.1. Dust Properties 



3.2. Lyc Transfer Solution 



The general solution to equation (3) is 





where 




(5) 
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With our assumed density profile n{r), and defining riy = t{R), Fiy = F{R), equation (4) 
becomes 
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I{Tiy,R/rc) 



where 



I{Tiy,Z = R/Tc) 



1+z 



X 



X 



exp 



-Tly 



X 



-1.5 



1 - + 



dx 



(l + r(,/3)-i(l + 0.72/^)-i . 



(6) 

(7) 
(8) 



The approximation (8) has a maximum error of 8% near z = \ and Tiy = 3 for 2; e [.01, 10''] 
and Tly e [0, 10]. 

A fraction of the incident Lyc pliotons is absorbed by dust and by the recombining H 
and He ions in the evaporation flow, and the rate of arrival of photons at the ionization 
front is reduced to Fiy = Siy/inR^q, where the absorption factor, q, is given by 



q {qe 



-Tly 



1) = 0.1i;I{Tiy,R/r,) 



(9) 



which derives from multiplying equation (6), by Siy/[4:7rR'^F^y]. Integrating (5) with (2) 
yields 



Tly 



criyC2 ^ 
1.5 X 1.6 aB q 



545 q V 10-21 cm2 



(l + i?/r,)i-5 
1 



(l + i?/r,)i-5 



(10) 



which is the dust absorption optical depth for Lyc photons through the freely expanding 
evaporation flow between the star and the ionization front. We may treat R/tc and Tiy 
as the two parameters defining the problem, and find q{Tiy, R/rc) and tp{Tiy, R/rc) from 
equations (9) and (10). With aiy = 10~^^a_2i cuP we obtain 

1 (1 + 0m8a_,jTiy) for R » r,, 
\ 0.61 (rjR) for R < r^, 



2.99 X 10^ Tfy e'ly 



1 + %/3 

54.5 T{Rlr^) 



1 + 



1 + Tlyl'i 



(11) 
(12) 



where 



T{RIt,) = {l + 0.72rjR)-'[l-{l + R/r,)-'-'Y' = 



for R^ Tf. 



(13) 



.93 for R < r-c. 

The photoevaporation flow is optically thick to Lyc radiation for g ^ 1, i.e. ip ^ 10. 
Known evaporation flows have values of -0 up to ~ 10^ and Tiy ~ 1 for, e.g., the "partially 
ionized globules" (PIGs) in the Trapezium cluster. 
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3.3. Str0mgren Number and Shock Speed 

Our second dimensionless parameter, the "Str0mgren number" 

^ ^ ^ 5 23 X 10- 4^. fi5!£5L!)' , (14) 

is the ratio of the width {Siy/ iTvR^)/ {aBn^) of an ionized Str0mgren layer of density no (in 
which recombinations balance ionizations), and the cloud size, 2rc- A cloud of initial density 
no and ^ > 1 (and 5tp > 1) that is first exposed to the ionizing radiation gets completely 
ionized by a rapid R-type ionization front that never stalls to develop a shock. When S is 
small however, only a thin surface layer is ionized before a shock forms and the shielding 
photoevaporation fiow is established. For large tp, the velocity of the shock is derived from 
the momentum jump condition. If we assume a strong shock, Vs ^ Vao, where Vao is the 
preshock Alfven speed, then pov'j ~ P2(^2 + ^i); using the identity Siy/4:TTR'^noC2 = {2ipSY^'^ 
and equation (9) with q ^ 1, we obtain 

Vs ~ C2 (99 Sy/^e-^'y/^ (1 + Tiy/3Y/\l + 0.72r,/i?)i/^ . (15) 

However, Vg cannot exceed C2, since otherwise the evaporating gas could not diverge as it 
leaves the IF and would accumulate between the cloud and the star, leading to increasing 
absorption and the deceleration of the shock (Bertoldi 1989b). 

If the shock compression is limited by the cloud's magnetic field, as is most likely the 
case for typical interstellar conditions, then the shock-compressed gas has a density 

rii ^ \pi no , (16) 

where vao = fyio,5km s^^ is the Alfven speed in the unperturbed, upstream gas; the Alfven 
speed in the compressed gas is vai — 2^^'^{vsVAoy^'^ ■ The clumps in molecular clouds show 
linewidths that suggest vao,5 ~ 1 — 3 (e.g. Bertoldi &: McKee 1992). The IF propagates into 
the shocked gas layer at the velocity 

viF = ^^^^^ ^ ^Ao 5"' e-^^yl' (1 + riyl^fl'il + Q.72rjRfl' , (17) 

following a derivation similar to that of (15). 

The pressure of the ionized gas that freely streams off the cloud surface would not 
suffice to drive a shock into the cloud if 5 ^ 10^^^;^o,5- However, such a situation would 
not occur in a cloud where the clumped, dense molecular gas was in pressure equilibrium 
with the tenuous interclump gas before an H II region developed. This can be seen from 
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the following consideration. When an H II region develops that engulfs a denser molecular 
clump, the pressure of the more tenuous interclump gas certainly rises as it becomes 
ionized. If the clump was initially in pressure equilibrium with the interclump gas, its 
enhanced pressure would now lead to compression of the clump. However, the pressure of 
a freely expanding evaporation flow would, if S ^ lO^'^w^Qg, be lower than the initial cloud 
pressure, and thus also lower than the pressure of the ionized interclump gas. We conclude 
that a steady, freely expanding evaporation flow could not form. Furthermore, it is easy to 
show that under such conditions the interclump gas is sufficiently dense to absorb all Lyc 
photons before they even reached the clump's surface, unless the star was located very close 
to the clump surface. Thus for very small values of S, the star's Str0mgren sphere in the 
interclump gas would be smaller than the clump. 

4. Merged Fronts: Ionization Erodes Dissociation Region 

Consider a molecular cloud of H nucleon density rii and effective radius of curvature 
rc that is suddenly exposed to a FUV (912-llloA) flux F^v = (S'„„/47ri?2)e-'^-" = xFii, 
where Fh = 1.21 x lO^cm^^s^^ is the flux for the diffuse UV spectrum of Habing (1968) 
in this wavelength interval, and = {(Juv/(^iy)Tiy is the UV optical depth due to dust in 
the ionized gas between the molecular cloud surface and the star; a^v is the effective dust 
attenuation cross section for FUV photons. For Ry = 5 the dust extinction cross section 
at lOOOA is (Text ~ 1-02 X 10 cm . We assume a dust albedo ^ 0.5 in the UV, and 
estimate that the effect of scattering is such that the effective attenuation cross section is 

— (Tabs + 0-5asca — O.TScTeaii !^ 0.76 X lO^^^cm^, and s = o^vloiy ~ 0.76/0.61 = 1.25. 
Note that s varies slightly with Ry: s ~ 1.20, 1.25, and 1.30 at Ry = 3.1, 5, and 6, 
respectively. The corresponding values for (cr„„/10~^^cm^) are 1.64, 0.61, and 0.38. 

When the dissociating radiation first strikes the cloud surface, a dissociation front 
propagates into the cloud at the velocity 

VnF = 2e U Fuv / m = 18.2 -L -A A km s-\ (18) 

0.5 0.15 rii 

where e ~ 0.5 is the fraction of the incident FUV photons which are absorbed by H2 , 
accounting for the fact that about half the photons are far from strong pumping lines 
and are thus absorbed by dust (Draine & Bertoldi 1996); fd — 0.15 is the fraction of 
H2 pumpings that lead to dissociation. As the DF propagates into the cloud, ever more 
UV photons are absorbed in the dissociated layer between the ionization and dissociation 
fronts (i.e., the region we term the PDR) by dust and the reformed H2. As a consequence, 
the DF slows down, and, if there were no erosion of the dissociated atomic layer by ionizing 
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radiation, would eventually stall when the PDR has attained an equilibrium between 
H2 destruction and formation. The initial propagation of a DF into a static cloud has 
recently been studied by Goldschmidt &: Sternberg (1995). 

Since the cloud is also exposed to ionizing radiation, an IF penetrates the cloud surface 
at the velocity vjp = Siy/Ai^R^qni. Initially q = 1 since no recombining dense gas shields 
the cloud surface, but after a short time ^ min(rc, -R)/2c2 ~ 4 x 10'' min ( r^c, i?pc) yr a 
dense evaporation flow is established and q is implicitly given by equations (11) and (12). 

Whenever ionizing and dissociating radiation strikes the surface of a molecular cloud, 
both the IF and DF begin to propagate into the cloud. When vip > vbf, the DF cannot 
separate from the IF and the ionization and dissociation fronts form a merged structure; 
a PDR will not be able to form. With the absorption factor q from equation (12), the 
condition vip > vdf becomes 



0.6 0.15 S„ ^ ^^^^ 



e fd S, 



uv 



^ _54.5 HR/rc) Tiy 



cr_2i l-\-Tiy/2> 



-(S-1)T(;, 



(19) 



For cr_2i = 0.61 and s = 1.25, the RHS of equation (19) (Fig. 2) attains a maximum value 
9.9 (9.2) for R/rc ^ 1(^ 1) at Tiy = 2.26. For the assumed values of fa and e, the IF will 
therefore always be faster than the DF when (Siy/Suv) > 9.9(9.2), but such high values of 
Siy/Suv are rarely the case for a stellar radiation field. When the LHS is smaller than 9.9 
(9.2), there will be a regime between a lower (provided the LHS exceeds 0.15) and upper 
Tiy between which the inequality (19) is satisfied so that the DF is faster than the IF and 
a PDR can form. Larger values of a_2i or s would enlarge the parameter space where the 
fronts merge (cf. Fig. 2). 

The low Tiy regime of overlap arises because with insignificant absorption in the 
evaporation flow, the IF will propagate faster than the DF when the flux of incident ionizing 
photons is higher than twice the flux of UV photons that lead to dissociation (estimated 
as ~ 7.5% of all incident UV photons). With increasing Tiy, the Lyc radiation is absorbed 
by both dust and recombined H - the UV radiation however is only absorbed by dust, 
resulting in an increase in the velocity of the DF relative to the IF. At high Tiy and s > 1, 
the difference between the effective cr„„ and uiy results in vip/vbf increasing exponentially 
with Tiy. The larger the value of s = Guv/criy, the lower the value of Tiy at which the IF is 
again faster than the DF (cf. Fig. 2). 

We computed the main sequence lifetime-averaged Lyc and FUV luminosities of O 
stars by combining the stellar evolution calculations of Schaller et al. (1991) with the 
non-LTE stellar atmosphere models of Kunze (1994). For stars of masses between 20 (09) 
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and 6OM0 (05), we find 

1.054V for .085 < 549 ^ 3 , (20) 

and this ratio levels off for higher masses to reach 2.2 for a star with initially 120AfQ (Kunze 
& Bertoldi, in preparation). For a molecular clump exposed to the ionizing and dissociating 
radiation of a single 05 star (for which 5*49 ~ 3), we expect the IF and DF to overlap for 
values of Tiy below 0.24 or above 4.2, whereas for a 2OM0 09 star with 5*49 ~ 0.085 overlap 
occurs only when Tiy < .018 or Tiy > 8.3. In Fig. 3a the vertical lines delineate the areas in 
parameter space to the very left and very right where overlap would occur. We are however 
not aware of any observed evaporation flows that are thick enough as to fall into the high 
Tiy regime of overlap. 

Note that for massive stars the surface temperature, and thereby the ratio of Lyc to UV 
emission can vary significantly during the evolution of the star. Young stars tend to have 
higher effective surface temperatures and lower surface metallicities in their photospheres, 
resulting in higher Lyc to UV emission ratios that gradually decline during the main 
sequence evolution. WR stars may also have higher than MS-average Lyc/UV emission 
ratios. 

If the radiation incident on the clump is due to a cluster of OB stars, the most 
massive star may dominate the Lyc irradiation, but the less massive stars may contribute 
significantly to the UV flux. An integration of the UV and Lyc emission over a stellar mass 
distribution similar to the field star IMF shows that the ratio of the integrated Lyc to 
UV emission still exceeds unity if the most massive star is earlier than 05. However, the 
relative strength of the UV and Lyc radiation on the cloud surface depends on the distances 
to the individual stars. 



5. Non-Equilibrium PDR: Transient Photochemistry 

A PDR is the neutral gas layer between the IF and DF where the abundance of H2 is 
low due to photodissociation by FUV photons. Stationary models of PDRs (which assume 
vdf = 0) postulate a steady-state balance between H2 destruction (mainly by UV-pumping) 
and H2 formation (mainly by grain catalysis) at each point. Grain catalysis of H2 is slow, 
with a timescale tgr = (27?3,.nn)"^ ~ 5 x 10^ yr/(nH/cm"^) (for Rgr = 3 x lO^^^cm^s"^; 
Jura [1975]). Hence the local rate per volume of UV pumping of H2 must everywhere be 
low. In stationary models of thick PDRs {x/n ^ 10~^ cm~^), the FUV is attenuated in a 
layer of dusty gas with very low H2 abundance before reaching the molecular region. 
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In an advancing PDR {vbf > 0), however, steady-state conditions do not apply, 
and there will be regions in which the local H2 pumping rate per volume is increased 
(Goldshmidt & Sternberg 1995). This increase is likely to heat the gas beyond its 
equilibrium temperature (London 1978) and can alter the chemical structure in the PDR. 

Let L be the width of the layer in which the transition from atomic to molecular 
H occurs, i.e., the width of the DF. The importance of "transient photochemistry" can 
be assessed by comparing the flow time through the DF, L/vjjp, with the timescale of 
H2 formation on grains, tgr- For this comparison to be meaningful we need to assume that 
the PDR exists, i.e., that the DF and IF are not merged and inequality (19) is not met. 
Then the IF and DF both propagate into the molecular cloud at the same speed, separated 
by a layer of mostly neutral and warm gas, the PDR. When the cloud is first exposed to the 
strong FUV and Lyc radiation, the DF may propagate more rapidly than the IF, but as the 
FUV flux arriving at the DF is attenuated by the dust and self-shielding gas between the 
IF and the DF, the DF slows until it reaches a velocity equal to that of the IF. If the DF 
initially propagates ahead of the IF-driven shock, the distance and column density between 
the IF and DF may in fact decrease in time (cf. §5.3; Hill & HoUenbach 1978). 



5.1. Width of a Stationary Dissociation Front 

In order to estimate the width of the DF, we first consider the structure of a stationary 
planar PDR that is illuminated by an FUV flux F^v = X^h in the wavelength interval 
912 — lllOA, where is measured at the IF, which constitutes the surface of the 
PDR. We define Fx = F^v/ {^X)uv as the average flux per unit wavelength in the interval 
(AA)„„ = 1110 - 912A= 198A. 

Let N2{x) be the column density of H2 from the surface to some distance x into the 
PDR, let nn = 2n(H2) + n(H) be the total hydrogen density, t{x) = auyN}i{x) the dust 
optical depth to FUV photons within the PDR, and let us define W = Wx/X = W^/u as 
the dimensionless total equivalent width of the H2 pumping lines. 

The dissociation rate at a column density N2 is 

CiN^) = h\Fx ^e-^ , (21) 
and a balance between H2 formation and destruction in a dusty PDR requires 

n(H2) C(A^2) = n(H) Rgr ■ (22) 

We have estimated the equivalent width for a gas with a moderate excitation temperature 
T lOOK and Doppler broadening of a few km/s, values typical in molecular clouds, and 
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found that for .V2 ^ 10^'' cm-^ dW/dN2 ~ 5.14 x lO^^^cm^ and ~ 0.14. Although an 
individual line's equivalent width grows oc nI'"^ for N2 ^ 10^^ cm , the total dissociation 
rate approximately follows a power law for 10^^ < N2 < 10^^ cm~^ (Draine & Bertoldi, in 
preparation): 

/ AT \ -3/4 

C(A^.) - 4.6 X 10-^1 X e-^ (y^h^) , (23) 

which in a simple form accounts for self-shielding, overlap (significant at N2 > lO^^'cm"^), 
and the variation of the dissociation fraction fd{N2) because of the lines' saturation 
at different column densities. The steady-state solution for the molecular fraction 
y = 2n(H2)/n(H) is now 



y 



2«»^"H , / ^] ( ,,,J,'; ) e*. . (24) 



C{N2) ■ Vxcm-3; V3 X lO-i^cm^s-V V10i4cm-2 



Let L = \d\iiy /dx\~li characterize the width of the H/H2 transition in a stationary DF. 
From equation (24) and \dN2/dx\y^i = riu/A, the column density nuL of the DF is then 
derived as 

nHL = a„j[l + r4/^]"' , (25) 

where 

= (/>o e-"^* , (26) 

^0 . 12.4 fiiiir^) (^^.^^r . (27) 



nn/cm ^/ \ R J \7.6 x 10 ^^cm^ 

and Tpiir = T{y = 1) is the UV opacity of the PDR up to the point where half the H is in 
molecular form. Numerical models of stationary PDRs (Draine & Bertoldi, in preparation) 
have 

Tj,dr ^ ln(l + 2.7</.o) . (28) 

Thus (f)o oc x/'^H is a dimensionless measure for the width of the PDR, which derives from 
a balance of H2 formation and dissociation. Now the column density of the DF is 



l + (2.7 + </.o-i) 



4/3' 



1 -1 



(29) 



Comparison with numerical models (Draine &: Bertoldi, in preparation) indicates that 
equation (29) tends to underestimate n^L, but is accurate to within a factor 2 for ^0 > 0.01. 
Note that this tendency to underestimate riuL for stationary fronts is compensated for by 
the fact that propagating DFs are in fact narrower than stationary DFs with the same (j)o. 
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5.2. Non-Equilibrium Condition 

Now we are able to assess the importance of transient photochemistry on the surface 
of a photoevaporating cloud where both the IF and DF propagate into the cloud at the 
velocity vip (eq.[17]). We now also take into account that the incident UV flux is attenuated 
in the ionized evaporation flow by t„„ ~ STiy before reaching the surface of the PDR. With 
equations (9), (15), and (16), and the identity {Siy/47rR'^noC2)'^ = 2tp5 we obtain 

^^^e-T,,(.-3/4) n for > r„ 

(1 + T;^/3)3/4 \ 0.85(re/i?)i/4 for R < r,. ^ ^ 

Advection will be important in the DF when VDF[dn{H2)/dx] > Rgrn}in{H) . Evaluated 
where y = 1, and using the estimate (29) for a stationary front, this becomes 

vnF > 4R,rmL = 1.6 || 7.6 x IQ-^^ cm^ ^ [i + (2.7 + <Ao-i)4/^] km s'^ (31) 

With Vbf = vjp = Siy/{A'wR?qni) [cf. eq.(17)], and q given by equation (12) and ni 
expressed through 00 ^ x/^i (6q-27), we derive that advection will be important when 

....o(l.53.,-)->0.100(^)-''(^) 

The LHS attains a minimum of 3.06 when 0o = 0.106, whereas for the typical dust properties 
we chose, the RHS attains a maximum value of (10.6, 9.8) (S'^^/S'/j^) (corresponding to the 
limits R/tc ^ I and <^ 1), when Tiy = 2.26. Thus the propagation of the IF and DF 
always has an important effect on the structure of the DF when the illuminating source 
has Sly/ Suv ^ 3.5. Since O stars tend to have flux ratios smaller than this value, whether 
or not propagation efi'ects are important depends on the speciflc values of 0o and riy. In 
Figures 3 we show for what range of values for </)o and Tiy condition (32) is fulfllled. Only 
inside the oval lines, which correspond to different values of Siy/Suv, is the DF slow enough 
that the chemical abundances reach near-equilibrium values - everywhere else transient 
photochemistry dominates the structure of the DF. In Figures 4 we display the same, 
but with (j)o replaced by the Str0mgren parameter S. In Figures 5 we display the same 
parameter space, but only for the cases Siy/Suv = 0.2, corresponding to a late O star, and 
Sly/ Suv = 2, corresponding to a very early or young O star, where in addition, we marked 
the regions of parameter values where the IF and DF overlap, where the cloud is instantly 
ionized by an R-type IF, and where initially the DF propagates ahead of the shock. 



6.62 (1^ 



1 + 



54.5:^ 



Tiy 



a_2i l + Tiy/3 
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5.3. Trapping of the DF in the Shocked Layer 

When an O star ignites in a clumpy GMC, the tenuous interclump gas is quickly 
ionized, leaving the dense molecular clumps exposed to the strong ionizing radiation that 
slowly evaporates and compresses them (Bertoldi 1989a, b; Bertoldi &: McKee 1990). The 
clumps may be located deep inside the H II region for a long time - they are not necessarily 
located at the edge of an ionization-bounded H II region - so that the attenuation of UV 
photons due to the tenuous H II region gas between the star and the clump may not be 
significant. The PDRs on the surfaces of such clumps differ therefore from those discussed 
by Hill & Hollenbach (1978), who investigated the interplay of a DF and IF at the edge of 
an ionization-bounded H II region. 

When a clump is first exposed to the strong dissociating and ionizing radiation, an IF 
rapidly ionizes a surface layer of width 25rc on the side of the cloud facing the star. Due to 
the high thermal pressure of this "Str0mgren layer," it expands both into the H II region 
and the cloud, thereby driving a shock into the cloud that assembles a layer of compressed 
gas between the ionization and shock fronts (Bertoldi 1989b). After a time ~ rc/2c2, a 
steady evaporation flow is established and the shocked layer's column density near the apex 
of the cloud surface grows steadily. The DF may initially advance ahead of the shock. This 
occurs when the shock speed is lower than the speed with which the DF can propagate 
into the unshocked cloud, i.e. when t'DF('^o) > ^'.s, which with equations (15) and (18) 
[substituting tiq for ni and making use of eqs.(9), (10), and (12)] requires 

The broken lines in Figures 5 delineate the range in parameter space where condition (33) is 
met. Even if the DF initially propagates ahead of the SF, it may not reach its equilibrium 
extent because the continual growth of the compressed layer and the enhanced H2 formation 
rate therein steadily decreases the equilibrium width of the PDR, until the shock engulfes 
the DF and conflnes it to the compressed layer between the IF and shock. While the DF 
propagates ahead of the shock, its speed relative to the unperturbed gas is about half the 
speed of the shock, and is therefore much higher than its speed after it is trapped in the 
shocked layer. The advective efl^ects discussed in §5.2 are therefore even more pronounced 
throughout the somewhat intricate initial phase. Once the PDR is trapped in the shocked 
layer, we found above that over a large region of parameter space, transient photochemistry 
due to strong advection may dominate the structure of the PDR. 
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6. Examples 

6.1. M17 Northern Bar 

To quote an example of a photoevaporating and dissociating globule where we would 
expect transient photochemistry to be important, we consider the globular structure 
observed in detail by Chrysostomou et al. (1993) and Gatley & Kaifu (1987: Fig. 7) in the 
northern bar of M17. Estimating a 25" = 0.27 pc radius of curvature and a projected 
distance of 1.9 pc from the group of one 04 and two 05 stars with total ^49 ~ 20 (Felli 
et al. 1984), we derive \ogip ~ 4.6 and log 5 ~ —2.3 — 2 log(no/10'' cm~^). Assuming the 
uncompressed density of the molecular globule hq 10'' cm^^ (the compressed gas appears 
to be at about 10^ cm^^), we find that the DF and IF may have separated (estimating 
Siy/Suv — 1-5 — 2), although their parameters fall well into the domain where transient 
photochemistry is important. 

An uncertaintly in this estimate derives from our neglect of absorption of UV and 
Lyc photons in the H II region between the globule and the stars. Because of the strong 
clumping of the ionized gas and the complicated geometry of the H II region (Felli et 
al. 1984) it is difficult to resolve this issue. However, neglecting intervening attenuation, 
we would expect the peak density of the ionized gas streaming off the globule's IF to be 
n2 = S/ {A.'wR^qv2) — 1700 cm^^, where we computed g ~ 76 and Tiy ~ 0.56, assuming 
dense cloud dust properties. Our predicted value of the peak density agrees well with the 
^ 1000 cm~^ estimate of Chrysostomou et al. (1992) from Br 7 observations, or with the 
density of 2400 cm~^ determined by Felli et al. (1984) for clump #11 based on VLA radio 
continuum observations. This indicates that the ionizing fiux arriving at the clump cannot 
be much lower than our estimate. 



6.2. Gum Globules 

The photoevaporating cometary globules in the Gum Nebula (cf. Bertoldi & McKee 
1990) are believed to be photoionized by either C Pup=HD66811 (04 I; Walborn 1972) or 
7^ Vel=HD68273 (WC8+09I; van der Hucht et al. 1981). Because of their large distance 
(~ 50 pc) from the ionizing star(s), the globules have small photoevaporation parameters, 
ip < 10^, and the outfiows will not be very dusty, with Tiy ^ 0.05. Thus the condition (19) 
for merging of the IF and DF will be satisfied for Siy/S^v ^ 0.5, which will certainly be 
satisfied for the ionizing star(s): these globules should therefore lack extended neutral H 
layers surrounding them. Note that even at d ^ 50pc CPnp or 7^Vel dominate the FUV 
starlight background. 
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6.3. Orion A: M42 

A very prominent PDR is the one illuminated by the Orion Trapezium stars at the 
surface of the molecular cloud from which the Trapezium cluster emerged. Here the ionizing 
stars are very close to the cloud surface, so that R <^rc and we need not worry about the 
actual cloud curvature, but only assume that the ionized gas is free to stream off the cloud 
surface. We assume our line-of-sight to be approximately normal to the cloud surface, which 
lies R ^ 0.2pc beyond 9^0vi C (Wen & O'Dell 1995). The ionizing flux from the Trapezium 
may actually have created a cup-like region (Pankonin, Walmsley & Harwit 1979; Wen & 
O'Dell 1995), so the planar boundary considered here is highly simplified. 

The most luminous star in the Trapezium, ^^Ori C, is difficult to classify: Walborn 
(1981) reported variations in the apparent spectral type from 06 to 04 during 7 consecutive 
nights! To estimate the total ionizing flux from the Trapezium, we note that the 23GHz 
map of Wilson & Pauls (1984) has a peak brightness temperature of 24K, very nearly 
centered on ^^Ori C; the brightness temperature falls to about 25% of the peak at an 
angular distance of 2.1'. A uniform density ionized sphere of radius Rs = 0.32pc (2.2' at 
the assumed distance of 500pc), density Ug = 2900 cm^^, and temperature T = 9000K, 
reproduces the peak surface brightness and produces 340 Jy at 23GHz, about 85% of the 
total measured flux of 400 Jy (Wilson &: Pauls 1984). The recombination rate in this 
ionized sphere would be 9.5 x 10^^ s~^, but the rate of emission of ionizing photons by 
the star must be substantially larger, since a large fraction of the ionizing photons will 
be absorbed by dust. The extinction curve toward ^^Ori C (Lee 1968) is characterized 
by Ry = 5.5 (Clayton, Cardelli & Mathis 1989), from which the extinction cross section 
at 912A would be estimated to be 9.65 x 10~^^cm^; taking the absorption cross section 
aiy ^ X 10~^^cm^, and Ue/nu = 1.08, the dust optical depth from center to edge is 1.33, 
so that absorption by dust is important. Following Petrosian, Silk & Field (1972), we 
estimate that ~64% of the ionizing photons are absorbed by dust; thus we would estimate 
Sly = 2.6 X 10^^ s^^. The model is obviously imperfect: M42 is not spherically symmetric, 
there must be a stellar-wind-supported central cavity, there is additional ionization beyond 
2.2', and the dust properties are uncertain. In the analysis below we take ^49 ~ 3 for the 
Trapezium stars, mainly due to ^^Ori C. 

If ^^Ori C radiates like a T = 45000 K blackbody with 5*49 = 2.6, then its luminosity 
L = 4.2 X IO^Lq, consistent with spectral type 05.5 ZAMS (Panagia 1973). The other 
Trapezium stars {9^ Ori A=B0.oV, 9^ Ori D=09.oV; Carruthers 1969) raise the total 
luminosity to L 5 x 1O^L0. We assume the Trapezium to lie a distance R ^ 0.2pc in front 
of the molecular cloud; if a fraction //r of the radiation from the star is converted to infrared 
emission, then the resulting integrated surface brightness of the cloud surface behind 
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the Trapezium would then be ~ Lfm/l&'K'^E? = 32//^ erg cm ^ sr \ in reasonable 
agreement with the measured ^ 10 ergcm^^ s^^ sr^^ (Werner et al. 1976; Drapatz et 
al. 1983). Note that if R = O.lpc, as favored by Tielens & Hollenbach 1985, L/IQ-k'^E? 
would exceed the measured value by an order of magnitude. 

From equation (1) we find ipR/rc ~ 6.7 x 10^. With aiy = 5 x 10"^^ cm^, equation (11) 
becomes ^ 

^ = 1.19x10^ ^'^^ ' x0.61{rjR) for i? < r-g , (34) 

1 + Tly/3 

from which we obtain Tiy ~ 1.65, q ~ 560, and n2V2 ~ 9.7 x 10^ cm~^ s~^. Taking V2 = I.6C2, 
we estimate the electron density at the ionization front to be ^2 = 5300 cm~^. The density 
in the PDR has been estimated to be rii = 2 x 10^ cm^^ (Tielens & Hollenbach 1985), 
and we assume an isothermal sound speed ci = 2.5 km s^\ corresponding to Ti ^ 1000 K. 
Rather than use the approximations (15) and (16), which are accurate only for a strong 
shock, we apply the momentum conservation jump condition across the ionization front: 

P2 [vl + 4 + {l/2)v%] = p, [vl + 4 + (1/2)^2 J ^ ^35^ 

where a transverse magnetic field is assumed. Since Vi = {P2/ Pi)v2, and = {P2/ Pi)vai, 
we may solve for vai- 

„2 _ {P2/P1){V2 + 4) - {P2/Plfvl - C\ 

1 - P2IPI 

The density ratio pi/p2 = 38 then implies vai = 3.4km s~^. Similarly, we may apply the 
momentum jump conditions to the shock front to obtain 

2 _ (pi/po)ci - eg + (pi/po - Pq/pi)v\J2 
1 - Po/pi 

If we assume a shock compression ratio pi/po = 2 (i.e., preshock density no = 1 x 10^ cm^^) 
and Co = 0.5kms^\ we obtain Vs = 6.5 km s^^ The preshock Alfven speed 
Vao = Vai{Po/ PiY^"^ ~ 2.4km s~^ is within the range 1-3 km s~^ believed to prevail in 
molecular clouds. The PDR gas would then be moving at a velocity (1— po/pi)^'s = 3.2 km s~^ 
relative to the preshock molecular cloud. This is in agreement with the observed velocity 
difference of 3 km s^^ between the C109q; line at vlsr = 11 km s^^ and the molecular gas 
at Vlsr ~ 8km s^^ (Jaffe & Pankonin 1978; Pankonin, Walmsley & Harwit 1979). 



Will the IF and DF be separate or merged? The initial speed ratio is 

viF ( Sly \ e''"y 



Vdf \SuvJ '^^fd.q 



(38) 
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With q = 560, r/y = 1.65, s ^ 1.3, and fa = 0.15, this becomes 
vif/vdf = 0.10(5'/y/S'„^)(0.5/e). If we assume that ^^Ori C radiates like a T 45000 K 
blackbody, then Siy/Suv ~ 2.4. This ratio is higher than that imphed by equation (20), but 
can be justified by the fact that very young stars tend to have higher effective temperatures 
and thereby higher Lyc to FUV flux ratios than averaged over their lifetime. Accounting for 
a population of later type O and B stars close to the cloud, we overall estimate Siy/Suv ~ 2. 
Therefore we estimate that vip/vdf ~ 0.2, a value which is insensitive to the chosen value 
of i?, but very sensitive to the ratio between the effective Lyc and FUV attenuation cross 
sections, s. With our chosen values, we find that the DF and IF will be able to separate 
until dust absorption in the PDR slows the DF to the speed of the IF. We also note that for 
an assumed distance R = 0.2pc from the star, the UV radiation field will be characterized 
by an intensity x = {Suv^^^^'^ / 4:'n'R^ Ffj) = 2.6 x 10'' - somewhat weaker than the UV 
intensity % = 10^ postulated in the PDR model of Tielens &: Hollenbach (1985). 

With equation (27) we estimate = 1-4. From equation (31) advection will be 
important for v^f > 0.30 km s~^. The velocity of the IF - and therefore of the DF - is 
Vi = 0.48 km s^^ (eq.[17] with ni = 2 x 10^ cm~^), implying that the DF is propagating 
rapidly enough that its structure is somewhat affected by transient photochemical effects. 
Note that the value we adopted for rii was based on modelling (Tielens & Hollenbach 1985) 
using a radiation field ~ 4 times stronger than our estimate. If rii is in fact smaller, vi 
would be larger, and the effects of advection on the abundances in the PDR enhanced. 
The strong absorption of Lyc radiation in the H II region and evaporation fiow makes the 
average energy of the Lyc photons arriving at the cloud surface very high, which would 
probably justify a lower value of aiy than the one we adopted. 



7. Discussion and Summary 

We have shown that existing models of stationary PDRs are often inapplicable to dense 
photodissociation regions around molecular clouds that are directly exposed to the ionizing 
and dissociating radiation from massive, early-type stars. Stars of early spectral types have 
sufliciently large values of Siy/Suv that, for a wide range of parameters, the DF is not able 
to separate from the IF. When this occurs, the IF and DF form a merged structure, and 
the UV pumping and photodissociation of H2 occurs to a considerable extent in gas which 
is also undergoing photoionization. Consequently, most dissociating FUV radiation as well 
as some of the Lyc radiation is absorbed by H2, not dust, and the efficiency for conversion 
of starlight into H2 line emission should be much enhanced over that of static, thick PDRs. 

Another consequence is that the UV pumping of H2 will take place in a region where 
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intense Lyman-ct radiation is present. If the Lyman-a is sufficiently intense, it may affect 
the IR spectrum of vibrationally-excited H2 , since the Lyman-a can selectively depopulate 
the ?; = 2, J = 5 and i) = 2, J = 6 levels (Sliull 1978; Black & van Dishoeck 1987). Whether the 
Lyman-o; intensities will be large enough to significantly affect the IR emission spectrum is 
unclear, but merits investigation. This may prove to be the explanation for the reported 
weakness of the 2-lS(3) line in the spectrum of NGC 6240 (Lester, Harvey & Carr 1988). 

112" will be present in the merged IF/DF, and it is possible that detectable levels of 
vibrational emission may be produced. The 112" will react with H2 to produce H3", which also 
may emit detectable levels of vibration-rotation emission. Theoretical modelling of merged 
ionization/dissociation fronts is currently underway (Bertoldi & Draine, in preparation). 

Even when the DF is able to propagate ahead of the IF, we find that the DP will 
often propagate too fast for stationary FDR models to provide a good approximation to 
the chemical and thermal structure of the PDR. When this occurs, we expect that the 
temperature of the molecular gas will be increased due to the increased rate of UV pumping 
followed by dissociation (about ~ 15% of the time) or collisional deexcitation (for densities 
n ^ lO^cm"^). If the density is high enough for collisional deexcitation to be effective, 
the DF may produce intense H2 vibrational line emission with a near-thermal spectrum. 
Stationary PDR models appear able to convert no more that 7 x 10^^ of the energy output 
from a star into emission in the H2 l-OS(l) line, in which case PDRs are unable to account 
for the intense emission observed from NGC 6240 (cf. Draine & Woods 1990). It appears 
possible that propagating DFs may attain much higher efficiencies for conversion of starlight 
into H2 line emission. 

We predict that in several famous dissociation regions such as M17N and Orion, 
non-equilibrium conditions persist in the dissociation fronts and their PDRs may not be well 
developed, implying smaller H I column densities than predicted by equilibrium models. 
Models of propagating ionization-dissociation fronts are required to theoretically describe 
these regions. Such models are currently under development and will be described in future 
papers. 

We thank D. Hollenbach for valuable discussions and criticism, and the referee, 
A. Sternberg, for helpful comments. This work was in part supported by the Deutsche 
Forschungsgemeinschaft and by NSF grants AST-9017082 and AST-9319283. 
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Figure 1 




Fig. 1. — Schematic rendering of conditions on the surface of a molecular cloud exposed to 
ionizing and dissociating radiation from a single star: the compressed gas layer between the 
ionization and shock fronts propagates into the cloud, driven by the pressure of the ionized 
gas that streams off the surface with an effective velocity V2 — 1.6c2, surface density n2 and 
sound speed C2. The dissociation front here is trapped in the compressed layer, but may 
initially propagate ahead of the shock. 
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Fig. 2. — Evaluating the condition under which the ionization and dissociation fronts merge: 
the right-hand-side of eq.(19) as a function of the Lyc optical depth of the photoevaporation 
flow, Tiy, for different combinations of aiy = cr_2il0~^^cm^ and s = auv/criy. For, e.g., 
Siy/Suv = 2 and dark cloud dust properties, the IF and DF merge for riy < 0.15 (cf. vertical 
line in Fig. 5b) or riy > 11. 
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Figure 3a 



photoevaporation parameter: log'^ 
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Fig. 3. — Cloud parameter space for R ^ Vc- (j)Q is a measure for the width of the PDR 
and strength of the UV radiation field, the photoevaporation parameter, -0, is a measure for 
the Lyc opacity, Tiy, of the evaporation flow. Equilibrium abundances are established in the 
PDR only inside the circular lines, which each corresponds to a different ratio of the Lyc to 
FUV flux incident on the cloud, (a) Dark cloud dust properties, cr_2i = 0.61, s = 1.25, (b) 
diffuse cloud dust, a_2i = 2, s = 1.20. 
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Figure 4a 
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Fig. 4. — Cloud parameter space for R ^ Tc, vaq = 2 km and (a) dark cloud dust 
properties, Ry = 5, (b) diffuse cloud dust properties, Ry = 3.1. The regime where 
the DF propagates too rapidly for equilibrium PDR models to apply (non-equilibrium) is 
separated from the regime where equilibrium abundances are established in the PDR by lines 
corresponding to different ratios of the Lyc to FUV flux incident on the cloud. 
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Figure 4b 
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Figure 5a p h otoevQ po ratio n parameter: log-^ 
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Fig. 5. — Cloud parameter space for R ^ Vc for dark cloud dust, a_2i = 0.61, s = 1.25, 
and vao,5 = 2: the Str0mgren number, S, is a measure of the dynamical impact of the 
photoevaporation on the cloud: clouds with S > 1 are rapidly ionized before a shock can 
form, shocks are driven into the cloud for smaller values of S, and for S ^ 10^^ the cloud 
is dynamically unaffected by the photoevaporation. (a) Siy/S^v = 0.2, corresponding to a 
~ 18Mq star, (b) Siy/S^v = 2, corresponding to a ~ 60Mq star. Only for cloud parameters 
within the "equilibrium" area are equilibrium abundances established in the PDRs, elsewhere 
transient photochemistry dominates. Clouds to the left of the vertical line in Fig. 5b have 
merged IF/DF structures [eq.(19)] and a PDR does not form. The region above the broken 
line is where the DF initially moves ahead of the shock front. For R <C Tc, one obtains 
almost identical diagrams if S is replaced by Src/R and -ip by ipR/vc- 
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Figure 5b 
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